Shear Viscosity of Hot Scalar Field Theory in the Real-Time Formalism 
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Within the closed time path formahsm a general nonperturbative expression is derived which 
resums through the Bethe-Salpter equation all leading order contributions to the shear viscosity 
in hot scalar field theory. Using a previously derived generalized fluctuation-dissipation theorem 
for nonlinear response functions in the real-time formalism, it is shown that the Bethe-Salpeter 
equation decouples in the so-called (r, a) basis. The general result is applied to scalar fleld theory 
with pure A(/>* and mixed g(f>'^+X(f>'^ interactions. In both cases our calculation conflrms the leading 
order expression for the shear viscosity previously obtained in the imaginary time formalism. 



I. INTRODUCTION 

The investigation of the transport properties of hot, 
weakly coupled relativistic plasmas is helpful for an un- 
derstanding of the collective effects in a Quark Gluon 
Plasma (QGP) produced in relativistic heavy ion col- 
lisions. The Kubo formulae of linear response theory 
101 provide a way to calculate transport coefficients us- 
ing Feynman diagram techniques in thermal field theory. 
The shear viscosity can be expressed as ||^,^ 
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To 



dt 

— QO 



dt'{n^-'{0)TT^,{^',t')y''' , (1) 



where tt^^ (/i, j^=0, 1, 2, 3) is the traceless viscous- 
pressure tensor which is imbedded in the energy- 
momentum tensor and in the comoving frame has only 
spatial components, TTy. For a scalar field tt^ = (SikSji — 
^S.ijSki)dk(t)di(j>,{i,j,k,l = 1,2,3). (...) stands for the 
thermal expectation value, and the superscript "ret" de- 
notes the retarded 2-point function which is related to the 
spectral function of the tttt correlator. Previous studies 
have shown that even at weak coupling a one-loop cal- 
culation of the spectral function of the composite 
field TTij is incomplete: the same finite thermal lifetime 
of the particles which is already required to regulate a 
pinch singularity in the one-loop result serves as a regu- 
lator for increasingly severe pinch singularities at higher 
order. Since the thermal width compensates for explicit 
coupling constant factors arising from the interaction ver- 
tices, an infinite number of multiloop diagrams is found 
to contribute at leading order to the transport coefficients 
. A nonperturbative calculation is required to resum 
this infinite set of leading order diagrams. 

In two pioneering papers, Jeon |^,|| derived within the 
Imaginary Time Formalism (ITF) a set of diagrammatic 
"cutting" rules, identified all diagrams which contribute 
to the viscosity at leading order, and then summed the 
geometric series of cut ladder diagrams to obtain the vis- 
cosity. This calculation was a remarkable tour de force 
which, to our knowledge, has never been repeated and 



checked. In this paper we will do so by reformulating the 
problem in real time using the Closed Time Path (CTP) 
formalism [p|-pl| . The present paper provides technical 
details for the short report given in ||] and generalizes 
that and other work to scalar fields with arbitrary 
interaction Lagrangians. In addition to giving an inde- 
pendent recalculation of the leading order shear viscosity 
in hot scalar field theory with 3- and 4-point interactions, 
the present work also provides a natural starting point 
for a future generalization to dynamical problems where 
the plasma is (slightly) out of thermal equilibrium. 

In the real-time formalism the resummation of an infi- 
nite number of leading-order contributions to the viscos- 
ity is done relatively easily by writing down and solving 
a Bethe-Salpeter (BS) integral equation for the 4-point 
function. With the help of the generalized Fluctuation- 
Dissipation Theorem (FDT) for nonlinear response fimc- 
tions jl^ and further simplifications arising from the 
so-called "pinch limit" (see discussion below), the BS 
equations for different thermal components of the 4-point 
Green function can be decoupled at leading order in the 
couphng constant [p|j7|p^. We use this here to derive a 
general expression for the nonperturbative calculation of 
the leading order contribution to the shear viscosity in 
scalar field theories with arbitrary interactions. Its ap- 
plication to scalar field theories with pure and mixed 
g(f>^+X(l)'^ interactions is then straightforward and repro- 
duces Jeon's ITF results in an economic way. A similar 
structure of the calculation is expected for gauge theories 
in the leading logarithmic approximation |14| . 

This paper is organized as follows: In Sec. p^I|we present 
the BS equation which resums all leading order diagrams 
to the shear viscosity, and we derive a general expres- 



sion for its nonperturbative calculation. In Sec. Ill 



we 



evaluate the kernel of the BS integral equation for pure 
A(/)^ and mixed gcfi^+Xcf)'^ scalar field theories. A short 
summary is given in Sec. ^S/\ 



1 



II. RESUMMATION OF LADDER DIAGRAMS 
AND A NONPERTURBATIVE EXPRESSION 
FOR THE SHEAR VISCOSITY 

A. CTP formalism in the (r, a) basis 

The Kubo formula (0) for the shear viscosity can be 
expressed in terms of the Fourier-transformed traceless 
stress tensor Wightman function |||,|^ , 

^ = li 0^™ n / '^'^ ^^"'^ ^""^^ ' 

where (3 is the inverse temperature. The Wightman func- 
tion (7rij(a;)7r*-' (0)) is identical with the (12)-component 
of the matrix propagator for the composite tt field in the 
CTP formalism. For any field "0 (composite or elemen- 
tary) the four components of this matrix propagator are 
defined as 



(3) 



where Tp represents the time ordering operator along the 
closed time path (corresponding, respectively, to normal 
and antichronological time ordering of operators with 
time arguments on its upper and lower branch), and 
0-1, 0-2 G {1, 2} indicate on which of the two branches the 
ip fields are located. The scalar field propagator will be 
simply denoted by A: 

Aa,a2{xi,X2) = -i (Tp (xi )(/)q2 (X2 )]) . (4) 

Following we define 

0q(x) = 0i(x)-(/>2(a;), (l)r{x) = i(0i(x)-|-02(a;)) (5) 
and the 2-point Green function in the (r, a) basis 

AqiQ2(xi,X2) = -i2''^'-^'^ {Tp[(f>ai{xi)(l)a2{x2)]) ■ (6) 

Here ai,a2G {a,r}, and rir is the number of r indices 
among {ai,a2}. Substituting Eq. (|^) into (|^) it is easy 
to show the following relations for the 4> propagator: 

Arr{xi,X2) = Ai2(a;i,X2) + A2i(a;i,X2) 

= -«({0(xi),(/.(a:2)}), (7a) 
Ara(a;i,a;2) = Aii(a;i,X2) - Ai2(a;i,a;2) 

= -*0(.T?~X§)([0(xi),(/.(x2)]), (7b) 

Aarixi,X2) = Aii(xi,a::2) - A2i(xi,a;2) 

= -^e{x°-x'i)mx2),H^l)]), (7c) 
Aaa(a;i,a;2) = 0. (7d) 

Obviously Araixi, X2) and Aar{xi, X2) are the usual re- 
tarded and advanced linear response functions. In ther- 
mal equilibrium the correlation function Arr{xi, X2) and 
the linear response functions satisfy the Fluctuation-Dis- 
sipation Theorem [nsl 



in momentum space where n{k^) = 1 / (e^'^ —1) is the 
Bose distribution. Eqs. (0) are inverted by 

Aai a2 (*^l 5 -^2 ) Aq,^q,2 (-^l ; -^2)^0101 Qa2a2 ' ('^) 

where repeated indices are summed over and 

(10) 



Qal — —Qa2 — Qrl — Qr2 — 



are the four elements of the orthogonal Keldysh trans- 
formation for 2-point functions pO| . For example, one 
can express in momentum space the (12)-component of 
the 2-point function related to the Wightman function in 
Eq. I) as 



A 



[Y\k) = i(A(r'(fc) - KTHk) + Ai7)(fc)). (11) 



Using the Fluctuation-Dissipation Theorem and taking 
into account A*.^{k) = Aar{k), we see that the (12)-com- 
ponent of the 2-point function A''^^) for the composite 
field TTij is purely imaginary in momentum space. Eq. (|^) 
can thus be expressed as 

77=^ lim A\Y\p) = -^ lim lmA\V'\p). (12) 




p+q 



FIG. 1. The shear viscosity can be expressed as an inte- 
gral over the 4-point single particle Green function, with two 
pairs of external legs joined to a point, in the limit of zero 
external momentum p. 

Since each composite field Tr^j is made from two single- 
particle fields (p, the shear viscosity ry is related to the 
4-point Green function of the <f) field, as illustrated in 
Fig. 1^. Substituting tt^ into Eq. ( p^ one gets 

/? [ d'^k [ d'^q 



5pO,p^oJ {2^Y j (2^)4 
X lmGii22{-k, p+k, q,~p~q)J^^ {q,~p-q) , (13) 

where the function Jijij), q) joins two (p propagators to a 
point (see Fig. |^), 



(14) 



Arr{k) - (l+2n(fc"))(A™(fc) - Aar{k)) 



(8) 



and Gii22{—k,p+k,q,—p—q) in Eq.(p^ denotes the 
Fourier-transformed (1122)-component of the (j) 4-point 
Green function. Its 16 CTP components are defined as 
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(xi,2:2,a;3,a;4) (15) 

= (-i)^ {Tp [<f>a,{xi)(t)a2{x2)(f>aAx3)(l>aAxi)]) 

where Gaia2a3a4{^i^ ^21^3^X4) is the corresponding 4- 
point Green function in the (r, a) basis, defined as [ pj]jT2| 

X (Tp [(/)Qj(xi)0a2(a;2)(/)a3(a;3)(?!)a^(x4)]) . (16) 

Generally, in the (r, a) basis the n-point Green function 
with only a indices vanishes: Gaa...a = p^ . 



B. The need for ladder resummation 

In the one-loop approximation the 4-point function of 
Fig. |l| reduces to a product of two 2-point Green func- 
tions Ai2 which, in the limit of vanishing external mo- 
mentum, have opposite momenta. The (12)-component 
of the 2-point function can be expressed in terms of its 
spectral function p For noninteracting particles, the 
poles of the two spectral functions with opposite loop 
momenta pinch the real axis in the complex energy plane 
from above and below, resulting in a pinch singularity 
for the integral over the energy circulating in the loop 
[ p^ . This singularity can be cured by using resummed (j) 
propagators with a two-loop self-energy in the denomina- 
tor P p^ , p^ . This resummation generates scalar quasi- 
particles with a finite thermal (coUisional) width and life- 
time which shifts the pinching poles away from the 
real energy axis. 

As shown by Jeon |^,|^, however, this regularization 
of the energy loop integral is not enough to produce a 
reliable leading-order result for the shear viscosity. We 
rephrase his arguments, using weakly coupled the- 
ory as an example. Using the two-loop resummed prop- 
agator in the basic one-loop skeleton diagram for 77, the 
propagator pair sharing the same loop momentum con- 
tributes a factor 0(1/ A^) to the shear viscosity. (We call 
this the "nearly pinching poles" contribution.) Clearly 
such terms involving "nearly pinching poles" are dom- 
inant in the weak coupling approximation. For multi- 
loop ladder diagrams with parallel rungs formed by one- 
loop diagrams connecting the two "side rail" propaga- 
tors (see Fig. ||), each pair of side rails shares the same 
loop momentum as the external momentum p approaches 
zero; correspondingly, its frequency integral generates a 
"nearly pinching poles" contribution ~ 1/A^ which just 
compensates the factor A^ from the two extra vertices. 

All multi-loop ladder diagrams with parallel rungs thus 
contribute at the same order as the simple one-loop dia- 
gram. On the other hand, multi-loop ladders with cross- 
ing rungs have different momenta on some of their side 
rails; since the integrals over the corresponding loop ener- 
gies are free from "nearly pinching poles" contributions, 
these diagrams are genuinely suppressed by additional 



powers of A^. Similar arguments rule out leading^rder 
contributions from all other multi-loop diagrams At 
leading order we must therefore only resum the infinite 
ladders with parallel rungs. In the OTP formalism this 
is achieved by solving a Bethe-Salpeter equation for the 
4-point Green function. 



p+k 



p+q 
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p+k p+q 
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FIG. 2. In \4>*' theory the summation of an infinite num- 
ber of planar ladder diagrams contributing to the shear vis- 
cosity at leading order results in a Bethe-Salpeter equation 
for the 4-point Green function. 

In the Keldysh basis, the BS integral equation couples 
the component 6*1122 needed in Eq. ( |l3| ) to three other 
components of the 4-point function: G1222, G2122, and 
G2222- The key technical issue is therefore whether and 
how the BS equation can be decoupled. In the following 
we will show that in the weak coupling limit this is indeed 
possible, and that this is most conveniently achieved in 
the (r, a) basis: all components of the bare 4-point ver- 
tex involving an even number of a indices vanish (see 
Eq. (|47|)), and the propagator Aaa = (see Eq. (|7d|)). 
When decomposing G1122 into a linear combination of 
the sixteen (r, a) 4-point functions an important simpli- 
fication arises from the generalized FTD ||l^ (see also 
]l3[ and recent work by Guerin [^) which relates these 
functions among each other, reducing the number of in- 
dependent 4-point vertices to seven. We will see that, at 
leading order in the weak coupling approximation, only 
one of these seven independent functions contributes to 
the shear viscosity and its BS equation decouples. Again, 
the generalized FDT plays an important role in the de- 
coupling procedure. 



C. Contributions to 77 in the (r, a) basis 



Let us now follow the steps of expressing Ec 
(r, a) components. From the KMS condition pC 
rives in momentum space JlT|p^ 



G. 



aia2a3a4 



(^1,^2,^3,^4) 

"G'aja2a3a4 (^1, ^2, ^3, ^4) 

{i\ai=2} 



(|13|) in 

one de- 



(17) 
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where all momenta flow into the vertex such that 
^1+^2+^3+^4 = 0, the star denotes complex conjugation, 
and ai — 2,1 for ai = 1,2, respectively. We can thus reex- 
press Gii22{-k,p+k,q,-p-q) as 

G1122 {~k,p+k, q, -p-q) 

= ~nip'^)(Gu22 + G*22ii)i~k,p+k,q,-p-q). (18) 

We combine this with Eq. ( |l5| ) to decompose ImGii22 in 
the (r, a) basis: 

Im G'1122 (-fc, P+k, q, -p-q) 

= -n(/)Im (^Gii22 + G2211) {-k,p+k, q, -p-q) 



— 4^(P )Im aaar^G aara Garaa G^aaa^^Grrra 
+ Grrar^Grarr^Garrrj {^k,p+k, q, ^p—q) ■ (19) 

The first four terms on the right hand side involve the 
fully retarded Green functions defined in pl| , Graaa, 
Garaa, Gaara, and Gaaar- They are part of the set of 7 
independent 4-point functions chosen in Ref. p2| which 



additionally includes G, 



Gn 



arrai ^arar^ 



and Gn 



The last 



four terms in Eq. ( |19| ) are not part of this set, but can be 
expressed through these 7 functions via the generalized 
FDT |l|]: 



Grrra{kl,k2, k'i,ki) = 



— NiGarra — ^2^(24)^ G^^jj^ — iVsTV^-g^j^ G*^^, 



Grrar{kl,k2,k3, k^) 



Grarr{kl,k2, k^, ki) 



N2NiGraaa ~ N2NzN\^^GI,,, 



Garrr{kl,k2, k'i,ki) = 



(fcl, ^2,^3,^4): 

(n\^^ + iV^ArgJ))^:,,,, - N^N.N^^^^Gl,,, 

rra ^'L^arur ' (34^ ^aarr (^1 ; ^2 ; ^3 1 ^4) ; 

N^NiGraaa " (^[23) + ^4 ^(24)^ ) ^araa " ^2A^3^^(2^3) - N2N4,N^^^^ G*^^^^ 

-iVaiVj^^jG:,,, - NiN^^l^G:,^,. - N.Gaarr] ikl,k2,k3, ki) , 

(1 + N2N:i + N2Ni + N^N^^G;,,, - N^N^Garaa ~ N2NiGaara ~ N2N:iGaaar 
arra 

N3G 

arar 

N2G 

aarr 

{ki,k2,k3,kA) . 



(20a) 
(20b) 
(20c) 
(20d) 



Here ^ N (kf) ^ l+2n{k^) and TV = , with 

A:?+/c5+fc§+fc2 = 0. Substituting Eqs. @) into Eq. (|l|) 
we express ImGii22(— fc,P+fc, 9, — (/) as 

Im G1122 {-k,p+k, q, -p-q) 

— Im [wAGaarr + ^Garar + CGarra + T^Graaa (^-^) 
+£Garaa + TGaara + QGaaaT\ {-k,p+k, q, -p~q) , 

where the coefficients A,B, ■ ■ ■ ,Q involve combinations 
of thermal distribution functions. In the relevant limit 



■ they are given, up to terms of order I3p''\ by 



A. 

V: 



f3p°n{p°) 



:C.^^((iVM) + (A.,^-l)), 
Pp'^n(j>'^ 



.(^Nk{N'^-l) + N,iNi-l) 



(22a) 

(22b) 

(22c) 
(22d) 

(22e) 



^l^i^N,iN^~l)~N,iNi~l)). (22f) 

Here Nk = N{k°) = l+2n{k°) and similarly for Nq. In 
the zero external momentum limit the explicitly shown 



terms are finite since limpo^o /^p" n(p'^) ~ 1, whereas the 
dropped terms of order (3p^ and higher vanish. Since 
£ = in this limit, Garaa does not contribute to the shear 
viscosity. 

The n-point functions involving only elementary fields 
(j) are symmetric under particle exchange: 



G, , ,Qj . . .aj . . . ( 

G. . .a, . . .cti . . . ( 



h ■ 



h ■ 



, k.j 



(23) 



Moreover, the factor Jij{—k,k)J^^{q,—q) in Eq. ( p3[ ) is 
invariant under each of the following changes of inte- 
gration variables: k^q, k^—k and q^—q. Using 
Eq. (|2^), both Graaa{—k,k,q,—q) under the variable 
change k^q and Gaara{—k, k, q, — g) under the variable 
change q —q become equal to Gaaar{—k, k, q, —q). Per- 
forming the same variable changes on the corresponding 
prefactors V and we can substitute in Eq. (O) 



[VGraaa + ^Gaara + GGaaarK-k, k, q, -q) 
^ lNkiN^-1) Gaaari-k, k, q, ~q). 



(24) 



Using Eq. (|3|) together with N{-k°) = ~N{k°) this is 
seen to be an odd function of k which integrates to zero 
in Eq. (13). Consequently, none of the four fully retarded 
4--point functions contributes to the shear viscosity. 
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Similarly, under the variable change q —^ —q we have 
Garar{~k, fc, q, -q) Garra{-k, k, q, -q) such that 

[BGarar + CGarra] {~k, k,q, —q) 

^ i ((A^I-1) + {N^-l))Garrai-k, fc, q, -q) 

^ ^{Nl-l)Garrai-k, fc, q, ~q) , (25) 

where the second substitution involves the variable 
change k^ q. We will see below that in the weak cou- 
pling limit this term does not contribute to 77 at leading 
order either, leaving only the contribution from Gaarr- 



D. Bethe-Salpeter equation for Gai 



Referring to Fig. ^ above and Fig. |5| further below, the 
BS equation for the 4-point Green function Gaia2a3a4 
can be expressed in the (r, a) basis as 

^ Gaia2ct3a4{ k^k^q^ q^ 

[zA„,„3(-fc)] [zA„,„,(A;)] i2n)^S\k-q) + 



HA 



ai0i 



-k)] [i^a,-nik)] 



(2^ 



X [i^Gp^^^a^aA-l, 1,(1^-1)] 



-?^/3i 71/3474 ( k,k,l, I) 
(26) 



As always repeated indices are summed over, and 
^/3i7i/3474 is the kernel (an amputed IPI 4-point vertex 
function) of the integral equation. For the component 

Garra WC havC 

Garra{ k, fc, q, q) 

= -iAra{k)Ara{k){2^)H\k~q) 

dH 

KrrPi-fi ( k, k,l, /) 



-Araik)Arr{k) j j^27r)4 

^Gf3^y^ra{-1, 1, q, ^q) 
-A™(fc)A™(fc) J 
'xG/34^4rai-l,l,q,-q) ■ 



-^^ra/34 74 ( k,k,l, I) 



(27) 



Here we used Aaj^a2i~k) = Aa^aiik) and Aaa(fc) = 0. 
In the second term we can reexpress A^.^ in terms of 
Ara and A^^ using the FDT (||). Since Ara{k) and 
Aar{k) are retarded and advanced Green functions, re- 
spectively, the poles of Ara{k)Ara{k) and Aarik)Aarik) 
in the complex fc" plane lie on the same side of the real 
axis. Hence, in the limit A ^ 0, they do not generate a 
pinch singularity in the integrand of Eq. (|l^). The term 
^ Ara(fc)Aar(fc) = I Ara(fc)P, On the Other hand, does 
generate a pinch singularity as A— >0; in the weak cou- 
pling limit, this term will thus dominate in Eq. (|l3| ) over 
the two other ones. With this "nearly pinching poles" 
approximation we can rewrite Eq. (|2^ ) as 

(28) 



Garra{-k, k, q, -q) - iV(fc°) | A™ (fc) | ^ 

r dH 



(27r)' 



-Krrl3i-ii{~k, k, I, ~l) G/3^j^ra{—l, I, 9, —q) 



Using N{—k'^) = —N{k'^) and the analogue of the symme- 
try relations ( |2^ ) for the amputated IPI vertex [|2|, we 
see that ^{N^ — l)Garra{—k,k,q,—q) is an odd function 
of k and thus gives a vanishing leading order contribution 
to the shear viscosity in Eq. (13). This leaves only Gaarr- 



d'^k 



i2n) 



jn(fc")(l+n(fc°))/,,,„,(fc) 



luiGaarri~k,k,q,-q)I^"'{q). 



(29) 



Here lTr,im{k) = —Jim{—k,k) is the notation used in 



E. Decoupling the BS equation for Gaarr 

We now show that the BS equation for Gaarr decouples 
from the other components of Gaia2a3ai in the "nearly 
pinching poles" approximation. From Eq. ( p6|) we obtain 

Gaarri k, k, q, q) 

= -Ara{k)Aar{k) [l(2^)H^{k^q) 

dH 



-Krrl3i-ti{~k, k, I, —l)Gi3^j^rri — l, I, 9, ~q) 



(2^)4- 

Ara{k)Aar{k)[l{2Tr)^S^ik^q) 



dH 



+Krraa{-k,k,l,~l)Gaarri~lJ,q,-q) }• (30) 

For the second equality we used K„-rr — which is the 
analogue for the amputated IPI 4-point vertex of the pre- 
viously mentioned relation Gaaaa = |l2j . At this point 
Gaarr Still couples to G^arr and G arrr ■ For later conve- 
nience we introduce a function M which is obtained by 
truncating two of the external legs of the 4-point Green 
fmiction G: 

[»Aai/3i(-0] [«Aa2/32(0] ^Pil32a3aA-^J,q,-q) ■ (31) 

It follows that 

Gaarr{~l, I, q, ~q) 

^ -Ara{l)Aar{l)Mrrrr{-l, I, 9, "?) • (32) 

As above we now use the "nearly pinching poles" ap- 
proximation, ignoring terms involving Ara{l)Ara{l) and 
Aar(l)Aar{l) relative to terms involving Ara{l)Aar{l)- 
Using Eqs. ( [7^ ) and (^, we can then write 
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Grarr{~l, I, ~q) 

= -N{l°)Arail)Aaril)Mrrrr{-l, I, 9, -q) 

= N{l°)GaarA-l,l,q,-q), (33a) 

^ N{f>)Arail)Aaril)Mrrrri'l, I, Q, 'Q) 
^~N{}°)Gaarr{-l,l,q.~q)- (33b) 

These leading order relations decouple the BS equation 
( ^ ) for Gaarr] we find 

Gaarri-k, fc, q, ~q) 

= ~Ara{k)Aarik) [ii27r)^S\k-q) 

'^^^ IC{~k,k,l,~l)Gaarr{-lJ,q,-q)\ (34) 



D'rik) = 



(27r)4 
with the kernel 



lC{-k,k,l-l) = 

(^N {I ) {Krrra^ Krrar) + Kr^ra^ ( k, k, I, I) . (35) 



F. An integral equation for the viscosity 

For further simplification we introduce the 2-point 
spectral density 

pik) ^ l[Ara{k) ~ Aar{k)] . (36) 

Writing the full retarded propagator Ara{k) as 

1 



^™^''^ k^-ReY.{k) + ilmJ:{k) 

and using Aar{k) = A*^(fc) we deduce 

2IniE(fc) 
(fc2 -RcI](fc))V (iniE(fc))^ ' 



pik) 



Ara{k)Aar{k) = |A™(fc)|2 - ^^^^ 



(37) 



(38) 



(39) 



2ImE(A;) 

Substituting Eq. (||) into Eq. (|3|) we get 

Im Gaarri^k, k, q, -q) = 

-^^^^lTaMrrrri~k,k,q,-q), (40) 

and the decoupled BS equation (|3^) can be written as 

lmMrrrr{~k,k,q,-q) = {27:)^6'^{k^q) (41) 

^ /C( — fc, fc, Z, — /)Ini Afrrrr( — ^, ^ <Z, — <?)• 



(27r)4 2IniS(0 
Defining further 



(2^ 



4™(q) Ini AW(-fc, k, q, -g) (42) 



and using it in Eqs. ( P9| ) and (|4^) we finally arrive at 
r^^l-J^^ n(fc")(l+r.(fc°)) ^^^^ ^^tl',^-^^^ ■ (43) 



lOy (27r)'^ 



Inil](fc) 



Here • = /Tr.im-D^™, and D^^lk) satisfies the integral 
equation (suppressing the tensor indices on Ij^, Dj^) 



DAk)^lAk)- J^IC{-k,k,l,^l) 



2IniI](/) ' 



(44) 



These last two expressions were previously obtained in 
1^,^. The present derivation, however, does not use a 
specific form of the interaction Lagrangian; within the 
"nearly pinching poles" approximation, Eqs. ( ^ ) and 
( ^ ) form a general result for all scalar field theories, with 
different forms of the interaction Lagrangian resulting in 
different kernels IC{—k,k,l,—l). We expect a similar in- 
tegral equation to hold for the leading order viscosity in 
pure (quarkless) QCD, with modified expressions for 
and /C. In the following Section, we evaluate /C in the 
CTP formalism for the and g4>^+X(t)'^ theories and 
show that the results agree with those obtained by Jeon 
Q in the imaginary time formalism. 



III. THE SHEAR VISCOSITY KERNEL IN 
SCALAR FIELD THEORIES 



A. Massless Xcj>'^ theory 



The Lagrangian for massless cj)'^ theory reads 



(45) 



and we consider it in the weak coupling limit A ^ 1. For 
this theory the kernel of the integral equation (see Fig. |^) 
can be expressed as 

1 f d^s 

Kpi1lPili{~k,k,l,—l) — I 7:7;T4 ^7172 7374^^/31/32/33/34 



2 J (27r)4 

X ^72/32 (■s)^73/33('5 + l-k) 



(46) 



where the prefactor i is the symmetry factor associated 
with the bubble connecting the two lines, and the bare 
4-point vertex in the (r, a) basis is given by [12] 



A, 



(47) 



na being the number of a indices among (ai, a2, 0^3, on)- 
The last equation implies that all bare vertices must in- 
volve an odd number of a indices (in particular X^rrr = 0). 
Using Eq. (13) and Aaaik) = 0, Eq. (H) gives 
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A2 /• d^S 



8 7 (2^)4 

Arr(s)Aar(s + Z — fc) + Ara(s) A^r (s+^ — fc) 



Krrari k,k,l, I) — I j;: 7 

8 J (Ztt) 



(48a) 



(48b) 



Aar{s)Arr{s + l — k) + A^r (s) A^q (s+Z — fc) 

A2 /• d^. 



(27r) 



Arr(s) Arr(s + Z — fc) 



+A™(s)Aar(s+Z-fc) + Aar(s)A™(s+?-fc) . (48c) 

Noticing that the s°-integrals over terms of the type 
Ara(s) Ara(s+'— fc) and lS.ar{s)l^ar{s+l—k) vanish be- 
cause the poles of these two terms are both on the same 
side of the real axis in the complex plane, we can make 
the following replacement: 



d'^s r 



(2^)4 
d^s 



(2vr) 



A^a(s)Aar(s + ?-fc) + Aar (s) A™ (s + ? -fc) 

Ara{s)-Aar{s)\ [Aar{s + l~k)- Ka{s + l-k)\ 



■ p{s) p{s+l—k) 



(49) 



Inserting Eqs. (^) into Eq. (g5|) and using Eqs. (|), (|3^) 
and (|49|), we express the kernel of the BS equation as 



(50) 



introducing for later convenience the shorthand notation 
l + n{f) f d^si d*S2 



f l + n{n f 



+ nikO)J (2^)4 (2^)4 
X {l+n{-s°))p{-S2) {2Tif5\s^-S2+l-k) . (51) 



Substituting this into Eq.(44) we obtain 



A2 f dH l+n(/0) p{l)D^{V} 
-n(fcO) Iml](/) 



d'^s 
(2^ 



n(s°)p(s)(l+n(.s"+fcO-/0))p(s+fc-0 . (52) 



The shear viscosity is obtained by inserting the solution 
of this integral equation into Eq. (Hh)- This result coin- 
cides with that obtained by Jeon pfin the ITF. 



B. Massive gcj)^ + A^* theory 



1. The (f> self energy 

We now investigate the contributions arising from an 
additional cubic interaction. For massless scalar fields, it 



renders the zero temperature theory perturbatively un- 
stable, causing a nonzero vacuum expectation value for 
(j). To avoid this we introduce a nonzero mass m for the 
scalar field and study 



^0 = ^(<9p(5 



m 



2 ^ 3!^ 



A 

4!' 



(53) 



in the weak-coupling limit {g /mf' ^\<^1. m, g, and A 
denote the physical mass and coupling constants at zero 
temperature, and Cc denotes the counter terms required 
for renormalization. The one- and two-loop diagrams for 
the self-energy in this theory are shown in Figs. ^ and IJ. 






(b) 



FIG. 3. One-loop self-energy diagrams in g(f>'^+\4>* theory. 

Diagrams |^(a) and |^(b) are independent of the ex- 
ternal momentum which does not enter the loop; they 
generate a purely real contribution 0{XT^—g'^T^/m^) 
to the scalar self energy. For T ^m/^/X this is much 
smaller than m?, and the viscosity can thus be calcu- 
lated using the zero-temperature mass m in all propa- 
gators 1^. For T^m/\/X, the thermal self energy be- 
comes comparable to and must be resummed, yielding 
quasi-particles with a thermal mass of order m^j^ ~ AT^ 
Jp|, p7|jl8| . For T^m/\/X the self-energy from diagram 
g(b) is then ^ — g'^T'^ / m^^; due to the weak-couphng con- 
dition g^ ^ Xm? this is smaller than diagram ||(a) by a 



'm^/(AT^)^l. At such high tempera- 



tures, the effects from the 4-point interaction thus dom- 
inate those from the 3-point interaction, and the theory 
becomes effectively a massless A(/)4 theory as discussed 
in the previous subsection. We thus concentrate on the 
interesting temperature range T~to/VA where on the 
one hand mass resummation is required but on the other 
hand 4)^ interaction effects are not yet negligible. 

The lowest order contribution to the imaginary part 
of the scalar self energy arises from diagram |3|(c). It 
vanishes on-shell (an on-shell particle cannot decay into 
two on-shell particles with the same thermal mass) but is 
non-zero and C(g^) for off-shell momenta. Its real part is 
smaller than that of Figs. ^(a,b) and does not contribute 
to the leading-order thermal mass The imaginary 
contribution vanishes for the side rails of the multi-ladder 
diagrams in the "nearly pinching poles" approximation 
which forces the side rail momenta approximately on- 
shell, but it plays an important role for single-line rungs 
in the BS equation generated by the exchange of a single 
scalar field (see Fig. ||(b) below). 
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(a) (b) (c) (d) 



FIG. 4. Two-loop self-energy diagrams in gcjy^ ->r\(f)^ theory. 

The leading order contribution to the on-shell imagi- 
nary part of the self-energy arises from the two-loop di- 
agrams m Fig.|. As shown in detail in g, it is ©(A^T^) 
(this includes the contributions from diagrams ^(b-c) in 
the weak-coupling limit (g/m)^ ~ 1). It results in 
an C'(1/(A^T)) thermal lifetime for the quasi-particles 
P, p7p^ . Just as in the simpler (f)'^ theory, the frequency 
integral over pairs of propagators from the side rails of 
the multi-ladder diagram sharing the same loop momen- 
tum is proportional to the quasi-particlc lifetime. The 
corresponding factor 1/A^ ~ w? /{g'^X) ~ {m/gY compen- 
sates for the extra coupling constants from the vertices 
associated with the rung. Thus again an infinite num- 
ber of ladder diagrams contributes to the leading order 
shear viscosity; they will be resummed through the BS 
equation. 



2. The BS kernel for gij)^ + Xfjy- theory 

The various contributions to the kernel of the BS equa- 
tion for the leading order shear viscosity in g(f)^ + Xcf)'^ 
theory are illustrated in Fig. |[ Whereas for the side rails 
of the ladder full (i.e. at least two-loop resummed) prop- 
agators must be used, in order to regulate the pinch sin- 
gularity in the integral over the side rail momentum, it is 
sufhcient at leading order to use bare propagators in the 
contributions to the rungs of the ladder shown in Fig. |^ 

. We will nonetheless formally evaluate them with full 
propagators, in keeping with the spirit of our skeleton 
diagram expansion, even though we will not make use of 
the implied higher order corrections which, for a consis- 
tent treatment, would also require the resummation of 
non-ladder diagrams and ladders with crossed rungs. 

Diagram ||(a) corresponds to the pure c/)'^ theory dis- 
cussed in the previous subsection. For pure 4>^ theory, 
one might think that the simplest kernel for the lad- 
der diagrams should arise from the single straight line 
rung (one-particle exchange) (not shown in Fig. |5|) . Us- 
ing naive power counting, this straight-line kernel would 
contribute only two powers of g from the explicit vertices 
which cannot even compensate for the "nearly pinching 
poles" contribution 1/A^ ~ (m/g)^ from the frequency 
integral over the pair of side rail propagators with the 
same loop momentum. If this were correct, in the weak- 
coupling limit each further single-line rung on the ladder 
would make the shear viscosity more singular. Fortu- 
nately, this fear is unnecessary: due to the nearly pinch- 
ing poles from the pairs of rail propagators, the loop in- 
tegrals are dominated by almost on-shell rail momenta. 



For non-zero values of the overall loop momentum k in 
Eq. ( |l3| ) this forces all rung momenta off-shell. In the 
(r, a) basis, only vertices with an odd number of a in- 
dices are non-zero; together with Aaa{k) = this implies 
that the kernel from the straight line rung is proportional 
to the spectral density of the propagator (see Eq. ( [s^ ) 
below) . For off-shell momenta, this spectral density van- 
ishes at tree level, so the single straight line rung with a 
bare exchanged propagator does not contribute. 

k q k k q 



p+k P+q P+k P+k p+q 




(a) (b) (c) (d) (e) 



+ 




(f) (g) (h) (i) (j) 



FIG. 5. Graphic representation of the Bethe-Salpeter 
equation for the 4-point Green function in g(f>^+X(f>'^ theory. 
Compared to the pure A0* theory shown in Fig. |^, the kernel 
giving the leading order contribution to the shear viscosity 
(shown in the second and third line) is more elaborate. With 
(g/m)^ ~ A, all contributions to the kernel are of order A^. 

A non-zero contribution to the kernel from single par- 
ticle exchange thus requires at least 1-loop accuracy for 
the propagator of the single-line rung; this is shown in 
Fig. ||(b). As discussed in the preceding subsection, the 
loop in diagram IHc) contributes an imaginary part of or- 
der g^ to the single-particle self-energy and spectral den- 
sity which is non-zero for off-shell momenta. Together 
with the two factors g from the vertices on the side rail 
this renders the kernel from Fig. |^(b) 0(^{g/m)^^ and 
thus of the same order as all other diagrams in Fig. |^ 
(which involve the exchange of more than one particle). 
Note that, whereas the exchange of a single bare propa- 
gator contributes zero to the BS kernel, its once iterated 
version, shown in Fig. ra(c), is non-zero and again of order 
0((g/m)4) [|. 

3. The individual contributions to the kernel 

The kernel from diagram |^(a) is given by Eq. (pO|): 

IC(-){-k,k,l,-l) = ^ Jd[si,S2]. (54) 

The other contributions are calculated in a similar way. 
We show the calculational process in some detail for dia- 
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grams ^(b) and ||(c) , but only mention a few key steps 
and list the final results for diagrams ^(d-j). 

In the (r, a) basis the bare 3-point vertex is given by 



<?QiQ2a3 ~ 23/2 (j' ( 1) °) ' 



(55) 



where Ua is the number of a indices among (ai, q;2, as). 

Exphcitly this reads CJrrr = graa = gara = ffaar = 0, gaaa = 
grra grar garr ^/v^- 

For the single-particle exchange Fig. |^(b) the general 
expression for the kernel can be written as 



^ii/SiasP-A k,k,l, I) 

= -«5aiQ2a35/3i/32/33A02Q2(fc-Z) 



(56) 



where A is the full single-particle propagator. With 
Eq. ( j^Ej ) we get 



-^rrrai k,k,l, I) — i ^ Ara{k I), 



(57a) 
(57b) 
(57c) 



Inserting these into Eq. ( pSD and using Eqs. (g), (|3^) and 
( ^ ) we obtain 

/C(')(-fc, fc, I, ^l) = -^[N{f) + iV(fcO-;")] p(fc-/) 
= -g2[7V(/°)-fiV(fc°-?")]|A™(fc-/)|^ImS(fc-0. (58) 

As promised, the kernel for diagram ^(b) is proportional 
to the spectral density (or the imaginary part of the self- 
energy). For off-shell k—l the lowest order contribution 
to ImS is given by diagram ^(c): 



ImS(fc-/) - l(E™(fc-0 - ^ar{k-l) 

= 4ii{^rris)[Aaris + k-l)-Ara{s + k-l)] 



(2 

+ [Ara{s)-Aar{s)]Arr{s+k-l)^ 

Using Eqs. (||), (^ and the identity 



(59) 



(60) 



n{zi+Z2^ \-Zn) = 

n{zi)n{z2) ■ ■■n{z„) 

[H-n(zi)][l+n(z2)] • • • [l+n{zn)]-n{zi)n{z2) ■ ■ ■n{zn) 



we can insert Eq. (|59|) back into Eq. (^8|). With a little 
algebra we find that the result can be written in the form 

/C('')(-A:,fc,Z,-0- Y j d[s^,S2]\Ara{si-S2)W (61) 
where we used the (5-function in Eq. ( pl| ) to replace 

k — l = Sl — S2- 



For diagram ||(c) the kernel can be expressed as 

-^ii''/3i73i53 (~^' ^' ^' ^0 ^ .9aict2a35';3i/32/33 57l7273ff'5l<52i53 



(27r)4 (*)^')'2<52 {s+l~k) 

X A035i(fc-s)A^jQ3(fc-s) 



(62) 



With Eqs. (p35|), (7d), and (jg) and dropping the terms in- 
volving Ara(k—s)Ara(k—s) and Aor(fc— s)Aar(fc— s), we 
obtain 



|A™(fc-s)rAar(s+^-A:) 



rrra 

4 J (27r)4 

X A„(s) + A^(fc"-s°)(A™(s)-Aar(s)) 

^rrari~k, fc, Z, — /) 

X A„(.s)+Ar(fc°-s")(A,,(s)-A,,(s) 

^rraai^k, k, I, —I) 

\Ara[k—s)\'' Arr{sA-l — k) 



(63a) 



(63b) 



g 



^ d^S , . .,2 



4 J (27r)4l 



A, 



.(s)+7V(fc"-s°)(A™(s)-A,,(s) 



(63c) 



Inserting these expressions into Eq. ( pq ) and using 
Eqs. d), (H) and @j we find 

/C(^)(-fc, fc, Z, -/) = / y S2] \Ara{k + S2)\^ 

- Y / c?[si,S2] [|A,,(si-fc)|' + |A™(si+0|'] . (64) 

In the second equality we used the 5-function in Eq.([5l|) 
to change variables Si <-^- —S2 and replace S2+k — > Si+l. 

Repeating this procedure we derive also the remaining 
contributions to the BS kernel: 



K,^^\-k,k,l,~l 
= g^ j d[si, S2 

JC^''\~k,k,l,~l 
= g'^Jd[si,S2 

= g'^ j d[si, S2 

IC^s\-k,k,l,-l 

K^^\-k,k,l,-l 
IC^'\-k,k,l,-l 



Re 



Re 



Re 



Ara{si+l)Aar{sl-k) 
Ara{si~k)Aar{si-S2) 



(65a) 
(65b) 
(65c) 



Ara{si+l)Aar{si-S2) 

^Xg^Jd[si,s2]ReAra{si-k), (65d) 

^ Xg'^ Jd[si,S2]ReAra{si+l) , (65e) 
= /C(j^(-fc,fc,/,-/) 
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2 



j d[si,S2] ReAra(si-S2) . 



(65f) 



In key to obtaining these results is, beyond using the 
relations already mentioned above, a judicious redefi- 
nition of integration variables in order to bring them 
all into the same form. We also exploited the equality 
Aaian{—k) = /^a2ai(k)- Noting that diagrams 5(e), 5(g) 
and 0(1) can be obtained by turning diagrams 5(f), 5(h) 
and|^(j) upside down, we can use the following relations 
to simplify the calculation of the last three: The symme- 
try of the diagrams (e,f) gives 



Substituting this into Eq. 
-N{f) we obtain 

IC^^\-k,k,l,-l) = IC^'''>{k,^k,~l,l) . (67) 

JC''^\—k,k,l,—l) is thus obtained from IC^'^\—k,k,l,—l) 



(66a) 
(66b) 
(66c) 

H) and using N{-f) = 



by simply replacing k - 
find 

IC^''\-k,k,l,-l) 
JC^^\-k,k,l,-l) 



-k and l^—l. Similarly we 



/C(3)(fc,- 
/C«(fc,- 



-1,1), 
-1,1). 



(68a) 
(68b) 



The last equation implies the equality of the kernels from 
diagrams ||(i) and ^(j). 



4- The final result 

We can now add the contributions ( |5^ ) , (|6^) , (p4|) , and 
( |65| ) to the BS kernel for the shear viscosity in -I- Xcjr' 
theory: 



(69) 



A + g2(^A™(si-fO + A™(si-fc) + Ara(si-S2)) 



After accounting for the different metric convention for 
the ITF propagator used in Ref. ^ which results in a 
relative minus sign to the propagator defined here in 
Eq. (37), this result agrees completely with the one ob- 
tained by Jeon in Eq. (4.43) of Ref. g. 

The integrand in Eq. ( |69| ) corresponds to the square of 
the tree- level two-body "scattering amplitude" in finite 
temperature g(l)'^ + Xip"^ theory. Starting from the result 
( ^ ) for pure (f)"^ theory, one can obtain the BS kernel 
for gcf)^ + theory by simply replacing the tree-level 
two-body scattering amplitude A in (j)'^ theory by the cor- 
responding more elaborate amplitude in g(l)^+X(f>^ theory. 
Substituting Eq. ( |69| ) into the BS integral equation (^) 
and solving the latter (which for pure 0** theory was done 
in 1^ and will not be repeated here), one arrives at the 
nonperturbative result for leading order shear viscosity. 
It is of order 77-rVA2-r3(m/c/)* 



IV. SUMMARY AND CONCLUSIONS 



For the shear viscosity in hot plasmas a nonperturba- 
tive calculation is needed, resumming an infinite number 
of ladder diagrams which all contribute at leading order. 
Using a real-time approach based on the CTP formal- 
ism |^[l0[|, we formulated this task in terms of a BS in- 
tegral equation for the 4-point Green function. In this 
approach, the BS equation is a tensor equation which 
couples the thermal components of the 4-point function 
among each other. We showed that, at leading order, 
it can be decoupled by using the (r, a) basis of Chou 
et al. rather than the familiar Schwinger-Keldysh 
basis ||8|-|10||. In the (r, a) basis, the leading order con- 
tribution to the shear viscosity r] is given by the single 
Gaarr Component of the 4-point function. Decoupling of 
its BS equation makes abundant use of the Generalized 
Fluctuation-Dissipation Theorem for non-linear response 
functions p2| . 

The present derivation of the leading order shear vis- 
cosity is considerably more compact than the previous 
ITF calculation by Jeon [B. It differs from the calcula- 
tion of Carrington et al. Rjby following a strict organiza- 
tion in terms of skeleton diagrams with full propagators 
and bare vertices which we find conceptually more ap- 
pealing. This approach clarifies a disagreement between 
Refs. j|] and Ref. |Q as to the type of diagrams to be re- 
summed in a complete leading order calculation: we do 
not find a need for including any diagrams beyond those 
already considered by Jeon 1^ . Within our approach we 
first derived a general expression for the shear viscosity 
which describes the high temperature limit of all scalar 
field theories, and then evaluated the specific forms of 
the BS integral kernel for the Xcj)^ and g0^-|-A0^ inter- 
action Lagrangians. Our results fully reproduce those 
obtained by Jeon Q in the ITF, providing an important 
check and simplification of that hallmark calculation. For 
(g/m)^ ~ A <C 1, both theories give a leading order result 
for the shear viscosity of order rj^T^/X^. 

Due to the need for resummation, the diagrammatic 
approach to transport coefficients in hot field theories, 
based on Kubo formulae, has developed a reputation for 
being extremely demanding and cumbersome. This is 
particularly true for the ITF formalism with its need 
for analytical continuation and the corresponding com- 
plicated cutting rules As a result, an alternative 
approach based on the Boltzmann equation in kinetic 
theory has recently become more popular At this 

point, the only published calculation of the leading-order 
shear viscosity in hot gauge theories |Q is based on that 
approach. But these calculations are not trivial either, 
and they cannot fully replace the kind of intuitive in- 
sight into the underlying physical mechanisms which is 
provided by a diagrammatic analysis. A recent computa- 
tion PSI of the color conductivity in QCD, using the ITF 
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diagrammatic approach, sucessfuUy reproduced an ear- 
lier result based on an effective Boltzmann equation. 
A still unpublished calculation of the shear viscosity 
in hot QED and QCD, using the same diagrammatic ITF 
approach, differs slightly from the result in |Q, while a 
recent paper by Valle Basagoiti [|6) which claims agree- 
ment between the Kubo formula and Boltzmann equation 
approaches employs a resummed 3-point vertex which 
apparently does not satisfy the Ward identity ||2^ . (The 
violation of the Ward identity may not affect the leading 
logarithmic result, but this requires further study.) It is 
our hope that the more compact real-time formalism in 
the (r, a) basis developed here may help to resolve the 
existing discrepancies and lead to a revival of resummed 
perturbation theory as a tool for computing relativistic 
transport coefhcients. We also believe that it may serve 
as a suitable starting point for perturbative calculations 
of non-equilibrium transport phenomena, although away 
from thermal equilibrium the simplifications provided by 
the Fluctuation-Dissipation Theorem no longer work and 
hence additional complications must be expected. 
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